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Abstract. We present 3-D modelling of the distribution and kinematics of the neutral hydrogen in the spiral galaxy ESO 123- 
G23. The optical appearance of this galaxy is an almost perfectly edge-on disk, while the neutral hydrogen is found to extend 
vertically out to about 15 kpc on either side of the galactic plane. The HI layer and the major features of the HI data cube can 
be successfully explained by a model dominated by a strong (about 30°) line-of-sight warp. Other models were tried, including 
a flare model and a two-component model, but they clearly do not reproduce the data. This is the first unambiguous detection 
of a galactic warp that has the maximum deviation from the central plane almost along the line-of-sight. No evidence for the 
presence of any companion galaxy is found in the HI data cube. Line-of-sight warps in edge-on galaxies are probably frequent, 
but escape detection as they are too weak. Moreover they may easily be mistaken as flares or 'thick disks'. A 3-D modelling of 
the HI layer as the one presented here is needed in order to distinguish between these possibilities. 
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1. Introduction 

The outer parts of the disks of spiral galaxies are frequently 
warped (Bosma 1991 Garcfa-Ruiz et al. 2002 1. Warps are best 
seen with HI observations, since the neutral hydrogen extends 
out to much larger radii than the optical disk. HI warps obey 
some general rules, one of them being that they usually start 
around R25 (Briggs 1990 1. The measured angle between the 
inner plane and the outermost observed HI ring spans a wide 
range of values, from a few degrees to a few tens in some ex- 
ceptional cases. Of the various possible orientations that warps 
can have, one is considered as "unfavourable", as the warp be- 
comes less visible: this is when the largest deviation from the 
central plane occurs along the line-of-sight. 

The origin of warps is still poorly understood. Different 
models have been proposed, for instance: the precession of a 
disk embedded in a flattened dark halo, misaligned with the 
disk itself (e.g. Toomre 1983 >, or the accretion of gas with an- 
gular momentum different from that of the disk (e.g. Ostriker 
& Binney 1989 1, but none is completely satisfactory. 

In an ongoing project to study the dark matter distribution 
of spiral galaxies a sample was culled from the list of Persic 
and Salucci ( 1995 1 selecting the best determined optical rota- 
tion curves; the selection criteria were: symmetry of the Ho- 
rotation curve, high angular extent, high HI flux and low lumi- 
nosity (Gentile et al., in prep.) The aim was to combine these 



optical data with HI measurements tracing the gravitational po- 
tentials much further out. One of these galaxies, the edge-on 
spiral ESO 123-G23, turned out to look remarkable in its HI 
distribution. Indeed the HI emission extends far away (up to 
about 15 kpc) from either side of the plane. Such a pattern can 
be produced by different phenomena including a flaring of the 
outer disk, a thick layer of gas or a warp along the line-of-sight. 
We have therefore modelled the HI layer of ESO 123-G23 con- 
sidering these possibilities. Our analysis shows that the correct 
explanation is a warping of the disk along the line-of-sight. The 
other possibilities are clearly ruled out. 



2. Observations 

ESO 123-G23 is an edge-on spiral galaxy located at a distance 
of 38 Mpc (Ho-15 km s -1 Mpc Not much was known 
about this galaxy except that it exhibits a slight warp in the 
outer optical disk (Reshetnikov & Combes 1999 1. We have ob- 
served ESO 123-G23 in the HI line for 2x12 hours using the 
Australia Telescope Compact Array (ATCA) 1 and combining 
the data from two array configurations (750D and 1 .5 A), pro- 
viding a beam of 30'.'6 x 25'.'0 and a velocity resolution (af- 
ter Hanning smoothing) of 6.6 km s _1 . In Fig.[2(upper panel) 
we display an optical - HI overlay showing the peculiar mor- 
phology of this galaxy: while the optical image resembles an 
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Fig. 1. Optical DSS image (grey scale) overlayed with the HI 
total intensity map of the observed data cube (top) and of the 
line-of-sight warp model (bottom). Contours are 1 x 10 2() cm" 2 
(pseudo-3<x defined similarly to Verheijen & Sancisi 200 f), 4, 
7, 10, ... x 10 20 cm 2 . The HI beam is shown in the top left 
corner. 

almost perfectly edge-on disk, the HI is much more extended 
perpendicularly to the major axis, displaying what looks like a 
very inclined central disk and diffuse emission at large (up to 
15 kpc) projected distances from the optical disk; this diffuse 
emission accounts for about 40 % of the total HI mass. The 
left column of Fig. |2] shows some channel maps of the HI data 
cube. Every fourth channel is shown, so the displayed maps 
are spaced by about 26 km s _1 . The central channel map, at 
a velocity of 2855.4 km s , is the one closer to the systemic 
velocity. The channel maps of ESO 123-G23 show the emis- 
sion from a central edge-on disk, plus fainter emission from 
regions above and below the plane. Such an emission is more 
visible in the central channels, close to systemic velocity. This 
may indicate either a lower rotation velocity of the gas above 
the plane or a lower inclination angle (i.e. towards face-on) of 
the outer parts of the disk. A warp in position angle is visible, 
especially at v = 2723.5 and 2987.3 km s , and it is quite sym- 
metric. ESO 123-G23 is reported to have no optical companion 
of size between R op t/2 and 2 R opt within a radius of 10 R opt 
(Karachentsev 1 1 99~3l . An inspection of the HI data cube did 
not disclose any companion in HI either. The total size of the 
data cube is 33' x 33' (dimensions of the primary beam) x 442 
km s . The position-velocity diagram along the major axis is 
shown in Fig.[3](upper panel), where the east side of the galaxy 
is slightly more extended but no clear kinematical lopsidedness 
is visible. On the p-v diagram is overlapped (dots) the rotation 
curve of the galaxy (see Section 3). 



We considered three phenomena that might mimic the ob- 
served HI structure and kinematics: a) a line-of-sight warp; b) a 
large flare in the outer part of the disk; and c) a two-component 
structure, thin + thick (slowly rotating) disk, similar to that 
proposed for NGC891 (Swaters et al. [T997b and NGC 2403 
(Fraternali et al. 2001 ). These phenomena have in common the 
fact of displaying emission away from the thin disk and are 
hard to disentangle except via detailed 3-D modelling of the HI 
layer. 

3. Models 

As mentioned above, we modelled the HI data cube in three 
different ways. All these models consist of a set of concentric 
circular rings (Rogstad et al. 119741 . each described by a set of 
geometrical and physical parameters, chosen such as to pro- 
vide the best agreement with the observed HI data cube. Some 
parameters were kept fixed for the three models. 

The parameters that the models had in common were: 
the radial surface density distribution, derived using the Lucy 
method (Warmels 1988); the HI velocity dispersion, chosen to 
go from 12 km s _1 in the inner parts to 7 km s in the outer 
parts (Kamphuis 1993 1, in order to have values as realistic as 
possible; a constant value of 10 km s would not have affected 
much the results, while values different from 10 km s _1 by more 
than 5 km s _1 would not give a good account of the observed 
data cube; the shape of the vertical distribution of the HI layer, 
chosen to be Gaussian with different scale-heights for the var- 
ious models (the choice of other distribution shapes does not 
significantly change the results); the position angle as a func- 
tion of radius, taken to slightly change counter-clockwise after 
^25 = 80'.'7 (~15 kpc), and then to significantly (10°) change 
clockwise further out (see third panel of Hig-EJ. By trying mod- 
els with different values of the position angle we estimate an 
accuracy of about 2° (errorbars in Fig. yjl, beyond which the 
disagreement with the observations is flagrant. 

The rotation curve was derived with the WAMET (WArped 
Modified Envelope Tracing) method (Gentile et al. 2002 
Vergani et al. 2002 1, a modification of the Envelope Tracing 
method (Sofue & Rubin 2001 1; a detailed description of this 
method will be given in a forthcoming paper (Gentile et al., 
in prep.). Only small corrections (< 5 km s~') were necessary 
to refine the rotation curve derived with our WAMET method. 
These corrections were made in order to give a better represen- 
tation of the observed HI data cube. This method, however, like 
the methods that derive the rotation curve from the position- 
velocity diagram, provides the product v rot sin(i). Therefore we 
have to account for the inclination, that can vary with radius 
within the same model. For the inner parts of the galaxy we 
used, in all the models, an inclination angle of 85°, a value sim- 
ilar to optical estimates (e.g. the Lyon-Meudon Extragalactic 
Database, LEDA, gives an inclination of 82.3°). 

In the warp model the inclination angle was lowered by up 
to ~ 30° beyond R25 in order to fit the observations (see second 
panel in HigEJ. With the same procedure we used for the po- 
sition angle, we estimated an error for the inclination angle of 
5°. The rotation curve was obtained with the WAMET method 
with the assumed values for the inclination angle. The FWHM 
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Fig. 2. Selected channel maps; every fourth channel map is displayed. Contours are -4.6, -2.3, 2.3 (2 cr), 4.6, 10, 15, 20, 30, 
50 mJy beam 1 . The columns represent the observed data cube, the line-of-sight warp, the flare and the two-component models. 
Heliocentric radial velocities are shown on right. 
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Fig. 3. Position-velocity diagram along the central major axis 
of the observed data cube (top) and of the line-of-sight warp 
model (bottom). Contours are -4.6, -2.3, 2.3 (2 a), 4.6, 10, 15, 
20, 30, 50 mJy beam -1 . The rotation curve computed separately 
for the two sides is projected onto the top diagram and the 
adopted rotation curve is projected onto the bottom diagram. 
The east and west sides are also indicated on the plot. 

of the Gaussian vertical distribution of the HI was chosen to be 
300 pc. 

The flare model differs from the warp model in having a 
constant inclination of 85° and a rotation curve obtained with 
such a constant value. The FWHM of the HI vertical distri- 
bution is 300 pc inside R25 increasing to a value of 25 kpc at 
R ~ 100". Such an unrealistic value is necessary to reproduce 
the HI distribution (lower contours) visible in the total HI map 
(Fig-E upper panel). 

Finally the two-component model is the sum of a thin and 
a thick disk in analogy to the previously cited models proposed 
for NGC 891 and NGC 2403. In this model the thin disk has a 
rotation velocity taken from the WAMET method with a con- 
stant inclination and a thickness of 300 pc (FWHM). For the 
thick disk (40% of the total HI mass) the rotation velocity was 
lowered by 30 km s _1 with respect to the thin disk and the thick- 
ness of the HI layer was taken to be 20 kpc (FWHM). Similarly 
to the flare model, this last unrealistic value is necessary to re- 
produce the total HI distribution. For both the flare and the two- 
component model a different choice of the inclination does not 
improve the agreement with the observed data cube. 

4. Discussion and conclusions 

In Fig.|2]we show some channel maps of the data cube com- 
pared with the channel maps of the line-of-sight warp model, 
the flare model, and the two-component model (from left to 
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Fig. 4. Illustration of the likely situation of the warp in 
ESO 123-G23 as seen by the observer. 



right). This figure illustrates why the flare model and the two- 
component model have to be rejected: apart from their some- 
what unphysical choice of parameters (although necessary to 
produce the vertical HI distribution), they cannot reproduce the 
observed channel maps. For instance, the flare model cannot 
reproduce at the same time the large north-south extension in 
the central channels and the relatively small north-south exten- 
sion visible in the extreme channels (v = 2723.5 and 2987.3 
km s _1 ). The two-component model is not able to reproduce 
the 'elbow' present at v = 2749.9 and 2802.6 km s _1 , which is 
also indicated on the receding side, at v = 2908.2 and 2960.9 
km s . 

The data are indeed fitted best by the model of a warp that 
shows maximum deviation from the plane close to the line-of- 
sight; Fig.|5]shows that such a model reproduces in detail all the 
key features of the observed HI data cube. See the upper panel 
of Fig.[5]for the adopted rotation curve for this model; the errors 
are the difference between the approaching and the receding 
side, considering also a minimum error of 2/sin(i) km s _1 . In 
the upper panel of Fig.[3]we show the derived rotation curve of 
ESO 123-G23 (separately for the two sides) projected onto the 
observed position-velocity diagram, and in the lower panel the 
adopted rotation curve (the average between the two sides) pro- 
jected onto the position-velocity diagram of the line-of-sight 
warp model. It is remarkable that the de-warped rotation curve 
of this galaxy is in good agreement with the Universal Rotation 
Curve of Persic et al ( 1996 1 while the raw curve, instead, had a 
very unusual decrease of 30 km s _1 from 7 kpc to 15 kpc, never 
found in galaxies with the maximum circular velocity of ESO 
123-G23. Fig. |4] shows the HI rings of the line-of-sight warp 
model as seen by the observer and Fig.[2(lower panel) shows 
the total HI map obtained with this model. The overall HI dis- 
tribution is well reproduced. We have also tried combinations 
of the various models, but the emission at projected distances 
of up to 15 kpc from the central disk can be explained only with 
the warp model, and combining the warp model with realistic 
values of an outer flare and a thick disk gives irrelevant changes 
to the model. 

In order to describe the geometry of the warp, as previ- 
ously done by Schwarz ( 1985), we define the angle 9 between 
the plane of the inner disk and the plane of the warped ring, 
and the angle ft between the line of nodes (at the intersection 
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Fig. 5. Parameters adopted for the line-of-sight warp model. 
From top to bottom, these are: rotation velocity, inclination and 
position angles, and the angles 8 and /3 (see text). 

of the plane of the ring and the plane of the inner disk) and the 
line of intersection between the plane of the sky and the plane 
of the inner disk. The uniqueness of this warp now shows up in 
Fig-Hi bottom panel: ft is always relatively close to zero, mean- 
ing that the maximum tilt of rings is almost exactly along the 
line-of-sight. We point out that a negative value of j3 means a 
position angle varying clockwise and that in unwarped sections 
of a galaxy f3 is of course not defined. The fourth panel from 
top of Fig. |3 shows the mono tonic ally increasing warp angle 9 
as a function of radius. 

There is, however, an ambiguity due to our assumption of 
a central inclination of 85° and the unknown true orientation 
of the disk: the inclination varies from 85° to about 55°, which 
could mean a variation of 30° or 40°, depending on the viewing 
angle. If we look at both the central disk and the outer rings 
"from above" or "from below", the warp angle is about 30°; 
if we look at one "from above" and the other "from below", 
the warp angle is about 40°. We plot the warp angles with the 
assumption of viewing the galaxy in such a way that the warp 
is as weak as possible (30°); considering the uncertainties in 
position angle and inclination, we estimate that the minimum 
warp angle necessary to reproduce the data cube is about 25°. 
Even with this assumption, this warp remains exceptional: the 
warp angle, measured at the last point, is very large. In fact 
it is larger for example than all the warp angles found in the 



systematic study of galaxy disks in HI by Garcfa-Ruiz et al. 
(2002). This is probably the reason why this warp is rendered 
visible: for smaller angles a warp oriented along the line-of- 
sight is almost invisible. 

We have thus found very strong evidence for ESO 123-G23 
to be an almost perfectly edge-on disk galaxy with a strong 
warp along the line-of-sight. This is demonstrated by the com- 
parison of a series of models with the observations. We there- 
fore report the first clear case of a galaxy with such a warp. 
We note that Becquaert & Combes (1997) interpreted the mor- 
phology and kinematics of the HI disk of NGC 891 in terms 
of a warp close to the line-of-sight. However, the fact that 
Swaters et al. (1997) interpreted their HI observations of NGC 
891 predominantly in terms of a thin + (slowly rotating) thick 
disk shows that the case was not as evident as the present one. 
Statistically, the frequency of occurrence of line-of-sight warps 
should be much higher, while most of them remain undisclosed 
owing to their small warp angles. The only way to discrimi- 
nate between different phenomena (warp, flare, etc.) is detailed 
modelling in the way described here. 
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